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We present detailed results from a recent microscopic theory of extreme correlated Fermi liquids, 
applied to the t-J model in two dimensions. We use two typical sets of band parameters relevant to 
the cuprate superconductors. The second order theory in the parameter A is valid in the overdoped 
regime where < n<0.7, where n is the particle density (related to hole density x by n = 1 — x). 
The calculation involves the self consistent solution of equations for an auxiliary Fermi liquid type 
Greens function and an adaptive spectral weight, or caparison factor, described in recent papers by 
Shastry (Ref. ([T]), Ref. ([5}). We present the numerical results at low as well as high T at various 
low densities in the normal phase. We display the (two) chemical potential(s), the momentum space 
occupation function m^, various energy dispersions locating the peaks of spectral functions, the 
optical conductivity, relaxation rates for quasiparticles, and the electronic spectral functions along 
various directions and with typical additional elastic scattering. The line shapes have characteristic 
shapes that are very similar to those from approximate treatments of the equations in recent work, 
and display frequently observed features such as the high energy kink, that depend on the bare band 
dispersion. 
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I. INTRODUCTION 



The t-J model describes the problem of extremely strongly interacting spin half Fermi systems, made difficult by 
the requirement of at most single occupancy of the lattice sites. It is the subject of many recent works in the context of 
the Cuprate superconductors and other correlated systems. This is also a very hard problem since standard methods 
such as perturbative expansions fail, thereby motivating a new strong coupling approach to the resulting in the theory 
of extremely correlated Fermi liquids Ref. (T), denoted by ECFL. Previous applications to the two dimensional cuprate 
systems of the ideas in Ref. (JJ have been approximate, but have given encouraging results. The spectral functions of 
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electrons that compare very well with the experimental Angle Resolved Photo Emission (ARPES) data as in Ref. ([2j), 
and lead to interesting predictions for the asymmetry of line shapes as in Ref. (3), as well as to unusual forms of 
electronic propagators in these systemsRef. (4). 

The recent work Ref. (5) provides a detailed set of equations following from the above ECFL theory, which are 
the main focus of this work. We will refer to this theory as (I) here and prefix equations of that paper with (I). 
The equations are to second order in the parameter A of that theory, related to the density of double occupation. 
Broadly speaking, these equations generalize the fully self consistent skeleton diagram expansion of the canonical 
Fermi liquid theory, as described in texts Ref. (6), Ref. (7), Ref. to the case of extreme correlations, and take care 
to satisfy the important shift invarianc^of the t-J model to each order. In this paper, we present the first detailed 
calculations from ECFL by solving the above equations to evaluate some thermodynamical variables, the spectral 
functions, ARPES line shapes and optical conductivity of the t-J model. Our calculations are valid in the high 
hole doping limit, i.e. the overdoped regime which is also the region of low electron density^ so that the second order 
theory in A is valid. Future work will be aimed at higher order calculations in A which would enable us to address 
densities close to optimal doping. The results compared with other approximations as well as a few experiments. 

While analytical methods beyond crude mean field theories have been in short supply, there is a valuable body of 
numerical results for the t-J model from exact diagonalization Ref. Cg), high temperature series expansions'^, finite 
temperature Lanczos method4^^^^, and variational wave function^I^H^. Noteworthy are the results of Ref. ([12]) from 
Prelovsek and co workers, who handle the series expansion in inverse temperature in a stochastic fashion, thereby 
obtaining results down to fairly low temperatures. Owing to finite size effects and the inherent nature of the high 
T expansion, the results from this theory, although broadly comparable to ours, seem more grainy. The Hubbard 
model at large U tends to the t-J model so the large U studies of this model are of interest. Quantum Monte 
Carlo methods, despite the difficulties associated with the sign problem, yield some valuable insights into the spectral 
features such as kinks'^. The dynamical mean field theory for the Hubbard modeP^^ gives a numerically exact 
solution in high enough dimensions of the model, and although strong coupling, i.e. U > W is challenging, there is 
impressive progress overall. The recent paper Ref. ([2T]) at strong coupling i.e. large U/t obtains detailed spectral 
functions that are roughly comparable to what we find here for the t-J model. 

The work proceeds as follows. In Section (II) we present a summary of the equations solved here from Ref. (5). In 
Section (III) we discuss the computational strategy and explain the scheme, using the Fast Fourier Transform method 
(FFT), so that the spectral functions can be computed efficiently. Section (IV) presents the detailed results of the 
calculation. Section (V) contains a summary and some concluding comments. 



II. SUMMARY OF THE 0{X^) THEORY 
In this theory the physical Greens function Q is obtained from: 

g{k) = g{k) f,{k) = . ^ ^^^J (1) 

Here /i(/c), termed the caparison factor, plays the role of an adaptive spectral weight while g{k) is a canonical Fermi 
liquid, derivable from a Luttinger Ward functional. In general both /i(/c) and g~^{k) can be expanded to a given 
order in the parameter A. Consequently we have two "self energies" denoted by $ and ^, which determine g{k) and 
/i(/c). The second order equations studied in detail in this work, are described by equations (1-83, 1-84, 1-85) for the 
Greens functions: 

9 

li{k) = 1-A-+A2 ^+A2^(A:), where {k) = {k,iujk), (2) 
^(^) = -^{^p^^k+q-p^^k^Sq^ Jk-p-uo) g{p) g{q) g{q ^ k - p) (3) 

p,q 

g-\k) = iuJn^^-£k-X^^{k) (4) 



£k 



3^2 \ 2 

l-An + A^^jsfc + A^] -Jk-q g{q) (5) 



£p -\- Sq -\- £k+q-p + Jk-p £p -\- £q -\- £k-\-q—p 
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This is written with shght change of notation ^{k) from (1-85), and we have introduced the effective band 

energy Sk in Eq. ([5| that gets a static contribution from shrinking of the bare energies Sk as weh as from the exchange 
energy J. Ah terms are understood to be correct upto O(A^), and possess corrections of O(A^) that are ignored here. 
We set the parameter A ^ 1 in the fohowing. 

The number of the physical electrons is fixed by the first sum rule: 

I = Y^Gik) e'-r^°\ (7) 

k 

while the auxiliary Fermions described by g are equal in number, and satisfy the sum rule: 

| = ^g(fc)e-"°" (8) 

k 

where = ^ ^ • We thus determine the two independent chemical potentials and uq, so as to satisfy the two 
equations (Eq. |7|) and Eq. (|8|) simultaneously. Before discussing the computation we note that dropping terms of 
O(A^) in Eq. ^ - Eq. ^ limits the regime of validity of these calculation to densities not too close to unity. From 
Eq. (p| we see that this theory would give a high frequency behavior of G ^ with cq = 1 — f + ^ , rather than the 
exact value cq = 1 — f . This slight error in the high frequency physics is a result of the expansion in A, however the 
low frequency physics epitomized by the Luttinger Ward sum rule is exactly obeyed to each order in A. Thus we have 
an error of ~ 25% in the high frequency spectral weight in this theory at n ~ .78, somewhat beyond where we can 
push this approximation. The O(A^) terms would extend the range of this approximation to higher particle densities. 
In the following, we will study the Greens functions and the optical conductivity as well as some thermodynamic 
variables for the 2-d square lattice t-J model. The detailed parameters are listed below in Subsectior jIV A[ 



III. COMPUTATION OF SPECTRAL FUNCTIONS 



A. Definitions 



Computationally it is expedient to employ a spectral function notation as described for example in Ref. (8). The 
Matsubara frequency object Q{k^iuJn) is analytically continued to the real axis and its spectral density written as: 

pg{k^ v) — — — Im [Q{k^ iou ^ v -\- if])] = A{k, v). (9) 

TT 

Here A(/c, v) is the same notation used in most experimental literature. The real part of the analytically continued 
function can be obtained by Hilbert transform 

Re g{k, v) = P.V. r ^^^^dy. (10) 

J-oo ^-y 

An analogous definition is given for spectral representation p^(k^v)^ p-^{k^u)^ p^{k^v) used for g, <l>, ^, etc, and 
hence the full set of equations above can be rewritten in terms of these spectral functions. Since ^ is a product as in 
Eq. ([l]) we note that within the O(A^) theory: 



pg{k,uj) =pg{k,uj) - 2 + T ^(^.^)J ^P^{k,uj) Re g{k,uj), (11) 

so the two sumrules Eq. (|8| and Eq. ^ can be written as: 

I = XI / dujp^{k,uj)f{uj) 

k 

^(l-f) = -E/ d^[p^{k,^) Re*(fc,w)+Reg(fc,w) (12) 
where fiuo) = (1 + exp(/3co'))~^ and /(cj) = 1 — fioo). The auxiliary spectral function is the usual Dysonian form 

Pg{k,uj) = 2 2- (13) 

{u ^ fi-Sk -Re^{k,uj)} + (ttp^) 
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The spectral functions for ^ and $ have the form 

p^{k,uj) = 1 dvidv2p^{p,vi)p^{q,V2)p^{p^q-k,iJi^V2-uj) X 

{/(^l)/(^2)/(z^l + - ^) + /(Z^l)/(Z^2)/(Z^1 + Z/2 - ^)} X 

+ Sk+q-p -^Sk^Sq^ Jk-p Uq) {Sk -\- €p -\- Sq -\- Sk+q-p + - {Jk-p + Jk-q) - } (14) 

p^{k,uj) = / djyidu2pg{p,ui)pg{q,U2)pg{p-\- q - k,ui -\- 1^2 - uj) X 

s pq J 

{/(z^l)/(^2)/(z^l + - ^) + /(Z^l)/(Z^2)/(Z^1 + i/2 - ^)} X 

(6p + Sk+q-p + + 5g + J/c_p — iXo) , (15) 

which differ only through the presence of momentum dependent factors. Here, Ng represents the number of sites on 
a finite lattice. We discretize the frequency integral into an fold sum, with a frequency step Acc; so as to cover an 
interval \uj\ < 15t. Typically we choose = 3000 such that Acc;=.01t. 

B. Flowchart of the Iterative process 

With these equations we can establish an iterative process which leads to a fully self consistent spectral function, 
pg. The self consistency loop proceeds as follows. 

1. Initialize all quantities to those of the Fermi gas: /jl = fiQ, p^ir = p^ = 0, uq is set to a plausible value of 2t. 

2. Build Pg from latest instance of /x, ^. 

3. Calculate from latest instance of pg and uq. Obtain the real part via Hilbert transform. 

4. Calculate new fi using a bisection root finder. 

5. Repeat steps 2-4 until /x and have converged to specified tolerance. 

6. Calculate from latest pg and uq. Obtain real part through Hilbert transform. 

7. Calculate J^^^ ^(^)g(^) and recalculate uq with some rootfinder. 

8. Return to Step 2 and repeat loop. Continue to the next step only when i^o has converged to specified tolerance. 

9. Calculate pg. 

The most computationally expensive step in this loop is the double integration for p^. If computed by a direct 
summation the computational time required will scale as N^N^. Furthermore, this slow step is on the innermost loop 
so it is repeated many times to find selfconsistent values of and uq. This leads to unacceptably slow convergence 
for any reasonable system size. Noting that the summation has the form of a convolution we can make use of FFT 
routines to calculate p^ with linear scaling in NgN^. This allows us to reach significantly larger systems and lower 
temperatures than would be possible by a direct approach. The next bottleneck in this flowchart is the calculation 
of the Hilbert transforms. These can also be made fast through a judicious use of FFT routines. Thus, by using this 
approach we obtain a scheme which can calculate the full frequency and momentum dependence of pg for lattices of 
substantial size, Ns ~ 2000 at temperatures as low as 30K. 

Is it useful to discuss the tolerances set on the Lagrange multipliers. The chemical potential is obtained to a 
precision of 10~^. This is significantly more accuracy than is required to satisfy the particle sumrule to within a tenth 
of a percent. However, we find empirically that this strict convergence criterion for /x can not be satisfied until the 
spectral function is also well converged. Thus, if fi has successfully converged to this tolerance we can be sure 

the is also well converged. The convergence criterion on uq requires the sumrule ^^(/c)^(/c) = ^(1 — f ) to 
be satisfied to less than 10~^. Again, this is overkill as it concerns the particle density alone. However, uq appears 
explicitly in ^ and ^. The chosen convergence criterion is such the final uq lands within .Olt of the exact value. This 
range is comparable to the smallest scales in our calculation, namely the frequency resolution Aoj and an implicit 
level broadening scale rj Acu. 
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To exactly calculate the spectral functions it is important to capture the entire range of the relevant frequencies. 
For non-interacting Fermions this requires a frequency window no larger than 8t. However, the spectral functions 
and have long tails which extend to much higher frequency even though the bare bandwidth is considerably 
decreased according to Eq. ([5|. Thus it is important to determine empirically what range of frequency is sufficient 
to capture the full support of the spectral function. As mentioned before, we employ a frequency grid which extends 
over the range |a;| < 15t, nearly four times the bare bandwidth, and find that this suffices to capture the support of 
all functions that arise. 



C. Fast Fourier Transforms for evaluating convolutions 

We next discuss the strategy used in the computation, where Fast Fourier Transforms (FFT) help reduce the time 
taken for the task. The equations for and have a form which is very similar to the particle-hole bubble diagrams 
familiar from a 2"^^ order perturbation treatment of the Hubbard model. 

E(fc)2n. + E(i,i)cxG(i,jfG(j,i) (16) 

pq 

where z = is a space time point. The convolution in Fourier space is a simple product in the space time domain 

and hence the real space version is advantageous. This is the well known core idea of the FFT technique, where the 
time savings arise since the Fourier transforms are performed in NlogN steps rather than N'^ (here N = Ns x N^). 
The ECFL "self energies", ^ and <l> have the same frequency convolution structure of T^{k)2nd which appears only 
through the frequency arguments of g. However, the and equations suffer from the presence of momentum- 
dependent decorations which render them not technically convolutions. Nonetheless we can use FFT routines to 
solve these summations. The strategy is to break up the integral into elementary pieces that do have the form of a 
convolution. We then avoid the need to do one large integral with quadratic complexity by doing many 70) small 
FFT's of linear complexity. 

To accomplish this we define several correlation functions which are similar to particle-hole bubbles. 

xoiQ) = ^g(9)g(9 + Q); xi(Q) = E^9s(«)s(« + Q) 

q q 

X2{Q) = ^^g+Qg(^)g(^ + Q); X3{Q) = Y.el^Qg{q)g{q^Q) 

q q 

(17) 

each of which is a full convolution in frequency and momentum and has a spectral function which can be calculated 
by FFT in linear time. 

With these correlation functions every ggg term (except one to be discussed later) found in Eq. ([3| and Eq. (|6| 
can be written in the form 

Sggg(fc) = F,{k)Y,F2{p)s{p)Xn{p - k)Fs{p - k). (18) 

p 

where ^1,^2, and F3 are each functions of momentum only and their arguments are carefully matched with the 
arguments of 5ggg, g, and Xn as they appear in the integral such that all factors fit the form of a convolution. In 
this way we can massage every term of into a convolution of one g and a Xn rather than three g's as originally 
written. There is one term in this problem which cannot be treated in this way because the argument matching of 



Eq. ( 18 ) cannot be achieved in such a simple way. This term looks like 

pq 

Nonetheless, this problematic term can be treated by the FFT approach if the factor Jq-k is broken up using an- 
gle addition identities. That the function Jk can be decomposed into a small number of Fourier components is a 
consequence of the locality of the exchange interaction. 

In defining the Fourier transforms, we need to extend the frequency functions to infinity, since it is only then 
that the frequency convolutions become products in the time domain. Recall that our frequency integrals have been 
discretized onto A^^; frequency bins which cover the support of our spectral functions. In extending the discretized 
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frequency summations to infinity, we follow the standard procedure of padding the A^^^ frequency bins with an equal 
number of frequency bins with value zero. By a simple exercise one can verify that padding finite data in this way 
allows an application of the periodic FFT in such a way that the result of the infinite transform is reproduced. This 
procedure is discussed in^^. No such considerations are required for the momentum sums which are by definition 
periodic and discrete, making them naturally suited to treatment by FFT. 

The Hilbert transform Eq. (10) is formally a convolution and can therefore be solved with the advantages of the 
FFT routines. Once again, however, we face the problem that this convolution is a non-periodic frequency integral. 
Furthermore, the Hilbert kernel ^, unlike other spectral functions with a compact support, falls off very slowly at 
large frequencies so the padding trick from the gx-convolution will not work well in this case. It is found that the use 
of FFT's to calculate a Hilbert transform will always introduce some error. Fortunately, this error can be controlled 
by increasing the length of padding used. In our code we use a frequency padding of SN^ for the Hilbert transforms. 
This relegates the error of the real parts of the various functions to very high frequency, far beyond the compact 
support of the spectral functions. The error introduced is therefore negligible. 



IV. RESULTS 



A. Physical Variables 

The program includes several parameters which can be varied. These include the tight binding bandstructure 
(through hopping parameters t, t' etc.), the spin coupling J, density, and temperature. For the parameters of the 
model, we study two important situations: 

• Case (A) the nearest neighbor hopping model with t ~ 3000K and J ~ 900K; 

• Case (B) with f = — At and the other parameters unchanged. 

Case (A) is the most natural one to study since it is minimal. However at the bare level, its Fermi surface near 
half filling remains closed around the F = (0, 0) point in the Brillouin zone. This is in contrast to angle resolved 
photoemission (ARPES) reconstructed Fermi surfaces in say BISSCO that open up, i.e. no longer surround the F 
point becoming hole like. For this reason we study the Case (B) as well, with a Fermi surface that closes around 
the r point for n< .55 and opens up at higher n. The Fermi momentum along the < 11 > direction is measured at 
n .85 in experiments on BISSCO as ~ .43A~^ and provides a useful check on the parameters. Within the Case 
(A) we find a considerably larger value ~ .55A~^ while Case (B) gives a better fit ~ .45A~^. While the Fermi volume 
is conserved by the theory reported here, the shape of the Fermi surface can, in principle, change since the self energy 
^{k) depends upon the wave vector ^, so these numbers may in principle change somewhat in the final spectrum. 
However, in this O(A^) calculation we find that the "bare" Fermi surface is preserved to high accuracy. 

We also note that for Case (B), the band density of states has a van Hove singularity (vHs) close to the bare 
Fermi energy for density n ~ .58. This ends up infiuencing the results for this case significantly in a range of densities 
around .6, and the contrast between this and Case (A) gives an idea of the role played by the van Hove points. 



B. Other parameters in the programs 



The code can furthermore be implemented on lattices of different size and spatial dimension. For a given choice of 
these parameters an appropriate choice must be made for computational grid. This includes the lattice size as well as 
the discretized frequency grid. We have looked at converged spectral functions for a wide variety of these parameters. 

The majority of the following results were performed on a square lattice which has dimension L x L where L=36 
and the lattice is taken to have periodic boundary conditions and translational symmetry. We therefore work in a 
momentum representation with an L x L sized k-grid of points kij = -^{hj) where 1 < i^j < L and a is the lattice 
parameter, taken to be a = 3. 82 A. The spectral functions have compact support which extends to \uj\<2Zt. We 
choose a frequency range |a;| < 15t = uJc/2 which is sufficient to capture the full range of the spectral functions. 
We discretize this frequency range in = 3000 bins each of width Auj = jf- = .Olt. Acj is the lowest resolvable 
frequency scale in the calculation so it is prudent to disallow any spectral features from becoming any sharper than 
this scale. Therefore we introduce the convergence factor rjmin = which will act as a lower limit on the width of 
spectral features. Thus in the Dysonian form of pg Eq. (13) we set P¥+ ^- -"-^ however reasonable to include 

larger r] in the calculation to imitate the effects of impurity scattering events in ARPES experiments. This parameter 
has been discussed previously in Ref. (2). 



7 




FIG. 1: Case (B) T = 130K. The density of dependence of the two Lagrange multiphers, /it and uq in eV. This density 
dependence is similar in character to that found for different choices of ( J, t^, T). The chemical potential approaches zero from 
below. This agrees with the expectation that /it = at half filling, uq is rising monotonically with density. The interpretation 
of this behavior is that with increased density, this O(A^) theory is unable to prevent double occupancy entirely. Thus the 
finite uq (in effect, a repulsive Hubbard U) must reenter the problem to penalize doubly occupied states. The scale of uq is 
comparable to the renormalized bandwidth leV). 



C. Thermodynamics 



In Fig. ([T]) we display the density dependence of the two Lagrange multipliers, fj, and uq^ namely the chemical 
potential, and the effective Coulomb energy for the auxiliary Fermions generated by the self consistency condition. 
The temperature dependence of /j. is shown below in Fig. ([2| revealing the temperature scale for the thermodynamics 
of the ECFL. For a Fermi liquid we expect that the chemical potential goes as /l^(T) = — bT'^ . We define a 
Temperature scale which sets the strength of the 0{T'^) term according to 




The temperature dependence of fj. is plotted in Fig. ([2| for Cases (A) and (B). Case (A) exhibits the expected 
quadratic dependence at low T at all densities up to the highest allowed by the O(A^) approximation. has been 
extracted and plotted in the inset as a function of the density n. However, in Case (B) the chemical potential sits 
very close to a van Hove singularity (vHs) in the very overdoped regime. Consequently, the quadratic form of /l^(T) 
is disturbed and we do not attempt to extract a for this case. The density of states of Case (B) is shown in the 
inset to illustrate the presence of the vHs at low frequency. Note the reduction of the scale of in Case (A) as we 
approach half filling, i.e. as the chemical potential approaches a vHs. 

The temperature dependence of uo{T) (shown as an inset of Fig. ([2|) differs from Ax(T) in several ways. Whereas 
the quadratic dependence of the chemical dependence occurs because of marginal changes in the occupation only at 
low frequency, uq is explicitly present in the form of ^ and ^ so small changes in uq are felt at all frequencies. In the 
high temperature limit uq acquires a linear temperature dependence (like fi) and the low temperature dependence is 
also linear, these are separated by a minimum at an intermediate T scale that it is larger in Case(B) than in Case 
(A). 



D. Frequency independent Variables 



We now proceed to study the Fermi surface in this theory through the momentum occupation function 

mk = ^g{k,icun)e'^-''^. (20) 

A drop in this function at the Fermi surface helps to locate the latter. This can be compared with the "Luttinger 
surface" defined by a sign change in Re^(/c,0) given in terms of the spectral function by 

Re gCk,0) = r (^^^^ (21) 
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FIG. 2: The chemical potential fj, versus T at n = .75 for Cases (A) and (B) (left and right respectively). An effective Fermi 
temperature T^, found from Eq. (19), is shown in Case (A) as a function of the density. In Case (B) the DOS of states is 
shown as an inset to illustrate the vHs. Note that the scale of variation of the chemical potential Afi ^ 15 meV in heating from 
T = lOOK to T = 300K, is quite large as compared to the changes that occur in weakly correlated matter, and are potentially 
observable in ARPES. A shrunken overall energy scale, seen clearly in Fig. ([5|, is ultimately responsible for this sensitivity. 
The temperature dependence of uq is shown in each case to be non-monotonic with a minimum at finite temperature. At high 
temperature uq rises linearly with T, similarly to /j.. 



At T = the Fermi surface in k space is traced out by Re Q~^{k,0) = 0, as dictated by the Luttinger Ward sum 
rule or the volume theorem. In Fig. ([3| is plotted the momentum distribution rrik at T = 130K for various densities 
of Cases (A) and (B) along three principle directions of the Brillioun zone. The Luttinger Ward zero crossings 
Re G~^{k^O) = are depicted by dashed vertical lines. There is a close correspondence between these crossings and 
the point where ruk = .5, similarly to that noted previously by Stephan and HorsclPin an exact diagonalization study. 
Since this correspondence is not on any rigorously firm basis, it is difficult to do more than to list the conditions for 
its approximate validity. Using high temperature expansions for the t-J model Singh and Glenistei^^ found the 
Fermi surface to be that of the Fermi gas by various criteria, and noted that the condition m/e^ ~ 0.5 is only satisfied 
approximately at high T. In Fig. Q we find that at higher temperature where the QP near the Fermi surface have 
been significantly broadened, and thus there are deviations from rrikp ~ 0.5 in Case (B), although Case(A) has much 
smaller deviations. In Fig. (|3| and Fig. Q, a point of considerable interest is the spill over of the occupation to the 
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FIG. 3: The quasiparticle occupation rrik is plotted along three principle lines of the Brillioun zone with T = 130K. The panel 
on the left shows Case (A) with a Fermi surface that is closed around k = (0, 0). Case (B) is shown in the right panel and has 
an open Fermi surface due to the presence of the finite t'. In each case the FS is the same as in the non-interacting problem. 
The Luttinger Ward crossing Re Q~^{k, 0) = is indicated for each density by the vertical dashed lines. For each density and 
each bandstructure the Luttinger Ward crossings correspond well with the condition rrik = 1/2. Case (A) shows much more 
variation as a function of density, because the gradient oi Sk is relatively small in the vicinity of the M-point. 



regions in k space that are unoccupied in the Fermi gas- as noted in various variational wave function studies of the 
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FIG. 4: The quasiparticle occupation function ruk is plotted along principle lines of the Brillioun zone for T = 605K. Unlike 
Fig. (Is]), the condition nikp — 1/2 is not satisfied closely in Case (B), although it is applicable to reasonable precision for Case 
(A). 



t-J model alread}ff^HI3. At each k value, since Eq. (20) can be written as a sum rule 



mk = J pg{k,uj)f{uj) duj, 



(22) 



we see that for momenta k > kp^ the spill over variable rrik provides an estimate of the spectral weight pg{k^uj) at 
energies below the chemical potential^ and should be of help in disentangling or resolving the "background" often seen 
but discarded in ARPES studies. 



E. Quasiparticle Dispersion and Zk 



We next study the location of the quasiparticle peaks. We can study three interesting functions to get a clear idea 
of the many body renormalizations that effect the dispersion relations in this theory. 



El 



/ 3n2\ 



rua 



efc -^ + Re <i>{k,Ek) 
max.[pg{k, oj) : w], 



(23) 



([5]) defines the "bare energy" with its static renormalization. Below we drop the J-term, which is quite 
0(1/20). Further locates the peak of the auxiliary Greens function g, and locates the peak of the 



where Eq 
small 

physical Greens function Q. This dispersion is complementary to E^ from Eq. (l23|. E^ in Eq. ( |23[ ) is the dispersion 
from the Energy Distribution Curve (EDC). Thus Ek and E^. give us the energies of the auxiliary Fermions and the 
physical Fermions. The theory has a strong asymmetry in the line shapes as noted in Ref. (3), and hence it is also 
interesting to track separately the locations of the left and right half maxima. In Fig. ([5| we illustrate the density 
dependence of the three dispersions in Eq. (23) for Cases (A) and (B). Furthermore, we provide insets which the 
bandwidths, W(n), of the various dispersions as a function of the density. Note the the bare dispersion of Ck is 2eV 
for both cases. In Fig. (|6| we plot the dispersions Ej. and E^ along with the asymmetric half maxima at a single 
density, n = .75. Near the Fermi surface we see that E^ ~ E^ but they differ far from the Fermi surface where the 
lineshapes become more extremely skewed. 

We next explore the difference between E^ and E^ near the bottom of the band (Fig. ([5|). The feature is reminiscent 
of the "high energy kink" , or more picturesquely "the waterfall" observed in the Cuprate compounds^^ and we next 
relate this to specific features of the spectral function. To do so it is useful to plot the Momentum Distribution Curve 
(MDC) dispersion EMDc{k) found by locating the maxima in k at a fixed energy uo 



k''{(jj) = max[p^(/c, cj) : /c], Emdc{^) = Inverse of k*{E). 



(24) 

The dispersion EMDc{k) can be compared with the energies in Eq. ([23|. In Fig. ^ the spectral function is depicted 
as a color density plot. Near the F-point where k = (0, 0) the QP becomes incoherent and the bulk of its spectral 
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FIG. 5: T — 130K. The vertical energy scale is in eV. The three dispersions defined in Eq. (23) are plotted along principle 
directions for three different densities. The insets show the bandwidth of the dispersions as a function of the density. The bare 
bandwidth is 2eV but each of these dispersions shrinks to a sma ller bandwidth. Here Sk is a k-independent scaling of the bare 
dispersion. The bandwidth renormalization due to Re$ in Eq. (23) is k-dependent, and so Ek has a different shape than e^. 
Note that Ek ^ El near the Fermi surface. However, in Case (A) El differs from Ek near the F-point for each of the densities. 
A lower inset on the left panel shows the evolution of real part of the denominator of g(/c, cj) with cj to illustrate the behavior 
oi El. A non-monotonic feature develops with a minimum near u — —.3eV which corresponds to the peak El. This peak 
disperses smoothly into Ek at higher k. However, Ek has a kink at small k where it disperses downward sharply. The pole Ek 
then acquires a wide peak which is short and broad. 




FIG. 6: Case(B) T = 130i^ We plot the physical quasiparticle dispersion El (Black solid) and the location of the two half 
maxima (Blue Dotted), as well as the auxiliary Fermion dispersion energy Ek from Eq. (23) (Red Dashed). The skew in the 
lineshapes of the physical particles leads to the nonsymmetric location of the black line relative to the blue dotted lines. 



weight is spread out to high negative frequencies. In this region the MDC peaks differ from and E^. This is the 
operational definition of the waterfall as seen in the Cuprates. The scales of the waterfall in this calculation are both 
slightly smaller than those seen in the Cuprates. The MDC curve recovers the scale of the bare dispersion e/.. 

Within this theory we can explore of various parameter dependences of the waterfall anomaly. In Fig. ([8| we show 
density plots of the spectral function with negative and positive t' /t. These correspond to Case (B) and a third case 
with t' /t = .4 which can be identified with the phenomenology of the electron doped Cuprates. In these plots we can 
see two major trends in the development of the high energy features of the spectral function. Firstly, the waterfall 
of the electron doped case is larger than that of Case(A) by a factor of 2. Case (B) has no measurable waterfall; 
the background at negative frequency is essentially featureless near the F-point. Meanwhile at positive frequencies, a 
similar feature develops near k = (tt, tt) as t' /t is decreased. Case (B) shows the best example of this inverted waterfall. 
Case (A) is an intermediate example where both features are present. The location, direction and magnitude of the 
waterfall is correlated with the points of high curvature of the dispersion. Regions where the QP band is relatively 
fiat show agreement between Emdc and Eedc- The waterfall is also visible, albeit with a decreased magnitude, at 
much lower density. 

We can also study the Fermi velocity along a given Fermi surface. In Fig. (|9| the angle dependence of the Fermi 
velocity is plotted over an eighth of the Brillioun zone for Cases (A) and (B). Case (A) shows a significant increase 
by a factor of 2 in velocity as we move from the antinodal to nodal directions. The absolute values of the velocities 
are roughly similar to those seen in experiments. 
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FIG. 7. L — 60, (n, T) = (.75, SOOiT). The spectral function of Case (A) is presented as a density plot to illustrate the structure 
of the background spectrum. We have overlaid the E^^ El, and EMDc{k) spectra as white, green and black points, respectively. 
In the low frequency region near /cf we see that the MDC, EDC, and Ek peaks coincide. Of particular interest however, is the 
region near k — (0, 0) where the MDC peaks occur at a significantly higher scale than Ek or E^.. Furthermore, in this region 
the EDC peak loses weight. Notice that the red peaks of the density plot, representing a sharp QP band, fade out to yellow 
near the F-point. This pronounced difference between the EDC and MDC peaks indicates that a waterfall- like phenomenon 
is present in this theory. There is also a new feature at positive frequency near k — (7r,7r) which looks waterfall- like. This is 
explored further below in Fig. ([8|. 



Finally we want to know the QP weight at the Fermi surface. In a conventional Fermi liquid the QP weight is 
defined by 



where S is the Dysonian Self energy and should be evaluated at {k^u) = (/cf,0). This definition does not 
immediately work for us as we do not have a conventional Dysonian form for our Q. To obtain an equivalent definition 
for Zk we define a new self energy according to the form 

- A (25) 



X - ^ X - S 

2 

where in this approximation A = 1 — | + ^ and x = icu -\- /jL — 6^. In analogy to the standard FL we now define 



as ' 



(26) 



Note that A differs from 1 — Rowing to the second order in A approximation made here. While it may be tempting 



to drop the term A from Eq. (26), it represents an important piece of physics in the larger context and therefore must 
be retained. To see this, note that the full spectral function of canonical electrons, e.g. in a large U Hubbard model, 
would have features at the scale of U that correspond to the upper Hubbard band, and are thrown out in the t-J 
model, where we have isolated the lower Hubbard band. Thus in a comprehensive canonical theory, the (low) value 
of Zk should be imagined to be compensated by a large piece I — Zk that is partly in the lower Hubbard band, and 
also partly in the upper Hubbard band (outside the domain of the t-J model). Numerics are performed at finite 
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FIG. 8: L = 60, (n, T) = (.75, 300ii^). The left and right plots show t' /t = .4 and t' /t = -.4, demonstrating that the background 
and the kink feature are sensitive to the bare band parameters. In these plots we overlay three dispersions to help with the 
visualization of the kinks. These are Emdc, Ek, and El (this is identical to Eedc) in respectively black, white, green. A 
positive t' ( usually associated with electron doping) results in the largest binding energy scale for the kink anomaly near the 
F-point. This features shrinks as t' is decreased. Therefore this calculation suggests that the relative magnitude of the waterfall 
in hole and electron doped Cuprates is a consequence of the sign of a result similar to the Hubbard model QMC results of 
Ref. (|18p . Simultaneously, a reciprocal kink at positive frequency, which we might call a waterspout, occurs around k — (7r,7r) 
with increasing magnitude for lower t'. This feature is therefore predicted by this calculation for future inverse photoemission 
studies. 



temperature and the T=0 value of is obtained by extrapolating the finite T data. It is found that falls quickly 
with density in both Cases (A) and (B) as plotted in Fig. (10). At the highest densities, is lower for Case (B) 
but it falls more abruptly at low densities in Case (A). It is clear that such small values of Zk speak of enormous 
backgrounds in the spectral functions, and below we capture the energy dependence of these backgrounds in detail. 



F. Self energies 



We now display the self energies that are involved in calculating the spectral functions. In Fig. (11) we show 
and p^. Both functions exhibit the uo'^ behavior very close to zero, as expected of a Fermi liquid self energy. 
Unlike conventional Fermi liquids (such as the order perturbation theory treatment of the Hubbard model), the 
magnitude of of the quadratic term is strongly k-dependent. From these functions and the associated real parts we 
can construct a Dyson self energy as defined earlier in Eq. (25). In Fig. (12) we plot the computed the imaginary 
part of the Dyson self energy p^,. It exhibits a similar magnitude and k-dependence at low frequency to that in 
p^. However, asymmetries begin to appear at intermediate frequencies. It is remarkable that at positive energy the 
function is considerably smaller than at negative energies, a feature that has already been noted in Ref. (3) and also 
in recent work in Ref. ([2T]) . 

The right panel of Fig. (12) shows low frequency fit parameters of p^] as a function of k. An appreciable low 



frequency asymmetry develops for unoccupied k. At higher frequency much larger asymmetries arise. Specifically, 
we see that for negative frequencies, a large scattering lobe extends outward to —uo ^ SeVwith considerable weight. 
This feature is larger for occupied wavevectors. No analogous feature appears for positive frequency. In Fig. (13) we 



show the T dependence of the single particle relaxation, obtained from the decay rate of the Dyson self energy at the 
quasiparticle pole ImS. The Case (B) is clearly unusual, but becomes understandable when we recall the proximity 
of the van Hove point near the Fermi energy for n ~ .58. Case (A) shows a crossover at about ~ 150 K to linear in T 
behavior, as compared to ~ 400 K discussed above in Fig. ([2|. These scales are lower than would be found in the 
2nd Qj.(^gj. treatment of the Hubbard model by roughly a factor of two (at this density). This difference arises because 
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FIG. 9: The Fermi surface and Fermi velocities for Cases (A) and (B) (left and right) with n = .75. The FS is marked out 
in a minimal octant of the square lattice. The gradients of the various dispersions {ek,£k,Ek) evaluated at the FS are called 
(vp^VF, vf)- In Ref. (23) the Fermi velocity of optimally doped BISSCO is measured as being nearly independent of the angle 
6 and has a value of 1.8 ± .2eVA, i.e. larger than ours (with our scale t ^ .25 eV) by factor of 3 to 4. 



1.0 




'*8.0 0.2 0.4 0.6 0.8 
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FIG. 10: The quasiparticle weight Zkp is plotted as a function of density for Cases (A) and (B). In each case Zkp vanishes 
with increased density, more quickly in Case (B), due to the vHs at n .58. The incoherent spectral contribution is therefore 
already ^ 4 times greater than the quasiparticle part at n ^ .7, and this ratio appears to increase further near half filling. 

of the bandwidth reduction which increases the important ratio T/W. Thus the rather early crossover is not caused 
by the large incoherent contributions to the Greens functions which arise in the O(A^). It is actually an effect of the 
0(A) terms. 

G. Spectral Functions and ARPES Lineshapes 

In this section, we present our results for the spectral functions within the second order theory to O(A^), argued 
above to be valid for low (high) enough particle (hole) density, say n < 0.7. In earlier work Ref. (2), we have compared 
the approximate results of the ECFL formalism with some phenomenological inputs, with the experimental data at 
somewhat higher particle densities n ^ 0.85. Ref. (2) found remarkably good fits with the line shapes in the nodal 
direction. We are content in this work to present the results from a microscopic calculation of ECFL with the simple 
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FIG. 11: Case (A) with {n,T) = (.75, ISOiC). The spectral functions for the two self energies $ and ^ , i.e. (left) and 
p^^f (right) in the two panels at several k points along the < 11 > direction vs. frequency. Both are roughly quadratic and 
symmetric at low frequency, but have a strongly k-dependent curvature. In the plot of p^, the minimum width 77 chops off the 
bottom of the low frequency minimum. 




FIG. 12: Case (A) with (n, T) = (.75, 130i^^). The Dyson self energy S can be inferred from the Greens function as in Eq. (25), 
and its spectral function displayed at several k points along the < 11 > direction vs. frequency. As with p^, has inherited 
a strong k dependence. We note here that the primary k-dependence at low frequency is seen in the softening of the quadratic 
coefficient with increased k. There is also an asymmetric term cubic term yos oc which grows in magnitude with increasing 
k beyond kp, as one sees from the panel at right, where we also see that in the occupied portion of the Brilliouin zone k < kp 
the curvature and cubic term have relatively little k-dependence. 



band dispersions in Case (A) and Case (B), at lower particle densities. We note below that the line shapes obtained 
here are of the same general nature as the ones in Ref. (2), although as one expects from a lower density situation, 
we find somewhat less dynamical asymmetry about zero energy. More detailed comparison with data near optimal 
doping with the microscopic ECFL theory must await the solution of the third or higher order equations where the 
criterion for validity discussed above (see para following Eq. ([8|) is satisfied more closely than here. 

Let us first note the local density of states at density (LDOS) n = 0.75 for both cases at low T in Fig. ([l4|. We note 
that the chemical potential is rather close to the peak in the LDOS for Case (B) at these fillings, and this proximity 
influences the results for this case significantly. 

We begin by displaying the spectral lines at very low density of particles in Fig. (15), where the spectral lines 
ioT k < kp begin to show a characteristic tails at high binding energy relative to the lines at k > kp^ which are 
more Lorentzian. Also interesting is the sharpening of the lines at higher momenta k ^ kp-, suggesting longer lived 
quasiparticles far outside the FS, in parallel to the results in Ref. (21). 

We next display in a set of figures Fig. (16) to Fig. (19) the spectral functions of the theory at low and high 
temperatures, at different densities, and along different directions. The first panel Fig. (16) shows the nodal spectral 
functions of Cases (A) and (B) at three different temperatures. The lines are quite sharp near kp but broaden out 
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FIG. 13: Case (B). The inverse lifetimes for QP's at the Fermi surface along < 11 > obtained from the imaginary part of the 
Dyson self energy, as a function of temperature, along with the momentum averaged inverse lifetime obtained from the optical 
conductivity as in in Eq. (28). The temperature dependence is more linear than would be expected of a Fermi liquid with 
an early onset, caused by the shrinking band width as seen directly in Fig. This situation in Case (B) appears even less 
FL-like due to the vHs. 




oj (eV) oj (eV) 

FIG. 14: For Cases (A) and (B) (left and right) with n = .75, we display the local density of states. The LDOS of the 
physical Q is shown in black while the bare DOS is the dashed red curve. The LDOS of the auxiliary g is shown as a dotted 
blue curve. The quasiparticle bandwidth is considerably more narrow and in each case there is a long tail extending to large 
negative frequency. Furthermore, it is clear that the tail of G decays more slowly (quickly) than g for large negative (positive) 
frequencies. This arises from the asymmetry produced by the caparison factor. 



rapidly away from kp- The insets give an idea of the change of spectral density with T. It is interesting that the 
detailed shapes at higher binding energy differ between the two cases, with Case (A) showing a more pronounced 
tendency for a rise with (binding) energy, giving rise to a minimum in some cases. This was also noted in Eq. ([5| in 
the discussion of the difference between Ek and . A feature of this kind has been observed in^ ^ and is reminiscent of 
asymmetrical features seen in ARPES experiment of LSCO. For Case (A) it is also noteworthy that lines for k > kp 
are sharper than the ones inside, and further, have tails that extend into < 0, thereby contributing to what is 
commonly relegated as "noise" or "background" in ARPES. 

Both Case (A) and Case (B) exhibit an interesting behavior in the QP scattering rates. This is seen in Fig. (17) 
over the first quadrant of the Brillioun zone. We plot the decay rate from the Dysonian self energy: F/e = p^{k^E^ 
as a function of k in Fig. (17). It is clear that in both cases there is a surface of minimum scattering rate which 
coincides with the Fermi surface. However, there is a second surface with k > kp^ where the scattering rates are a 
local minimum. This feature is particularly strong in Case (B) and can be seen clearly in the lineshapes of Fig. ([l6| 
and even at the very low density spectral function seen in Fig. ([l5|. In Case (A) the feature is less pronounced but 
is nonetheless present. 

We next discuss the role of elastic scattering in the lineshapes following Ref. (2). As discussed in^ the parameter 77 
can be used to model the effects of impurity scattering upon ARPES spectra and has invoked to distinguish between 
laser ARPES and synchrotron ARPES. Using fits from^ we arrive at 77 = .032eV and 77 = .12eV as being the typical 
values for describing the ARPES lineshapes of laser and synchrotron ARPES respectively. It is therefore interesting 
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FIG. 15: Case (B) with (n^T) — (.35, 130X). The physical electronic spectral function pg at several k points along the < 11 > 
direction is plotted, showing the evolution of lineshapes across the nodal Fermi surface. At this rather low density, the high 
frequency tails decay more rapidly than they do at higher densities. The inset shows the ECFL density of states overlayed 
with the Fermi Gas DOS for the same parameters. The difference between the Fermi Gas and ECFL at this density is much 
less than was seen in Fig. (14). Note also that the rather sharp QP's kp are visible at this low density, they persist to 

higher densities as seen below. 



to see how the line shapes discussed above in Fig. (16) evolve under the influence of the elastic scattering term. In 
Fig. (19) low temperature spectral functions are presented with the additional t]. The effect of adding a finite r] is 
similar to increasing the temperature, resulting in an increase of the linewidths. However, increasing temperature 
effects the linewidths more strongly near kp while rj broadens the lineshapes in a more broad and less selective fashion. 



H. Optical conductivity 

The optical conductivity, cf(Q)^ is computed within the bubble approximation here by throwing out vertex correc- 
tions: 

Recj{^) = ^Y.'^l f p^{k,uj)p^{k,n^uj){f{uj)-f{n^uj))duj (27) 

k 

The imaginary part of the conductivity can be obtained by a Hilbert transform of the real part. In this purely t— J 
calculation we must be careful how we combine the real and imaginary parts because the true imaginary part, as 
opposed to the one computed in the t-J model, will include contributions from the upper Hubbard band and from 
charge transfer processes which can be appreciable at high frequencies. For our current purposes we will simply discuss 
QP lifetime effects of two kinds. First we have I/to- which describes a characteristic lifetime which can be extracted 
from the low frequency curve of cr{Q). To do so we assume a roughly Lorentzian form and define l/r^r by 



i 



a{uj)/a{0)duj = 7r/4. (28) 



This is a momentum averaged quantity. Alternatively we can look at the momentum resolved scattering lifetimes 
defined as the inverse width of the ARPES lineshape at the Fermi momentum. These scattering rates are displayed 



in Fig. (20). The ARPES derived 1/r turns out to be essentially isotropic with differences of < 1% between nodal 



and antinodal Fermi surface so only the nodal r is displayed. 



In Case (A) the 1/r curves taken from ARPES and a{uj) in Fig. (13) rise quadratically at low temperature in 



accordance with the standard Fermi liquid picture. We note that the computed conductivity match up on an absolute 



scale quite well with experimental results in the overdoped regime. For example in Fig. (21) we display the Kea{oo) 
curves for Case (B) along with optical conductivity measurements published by Puchkov et al in^^ for an overdoped 
Thallium compound. 



V. CONCLUDING REMARKS 



In this work, we have presented a systematic low density expansion results for the t-J model, using the newly 
developed formalism of extremely correlated Fermi liquids^ discussed in Ref. (1) and Ref. (j^. The expansion in 
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FIG. 16: The physical electronic spectral function pg at several selected k points along the < 11 > direction are plotted, 
showing the evolution of lineshapes across the Fermi surface. Cases (A) and (B) on the left and right respectively with 
temperature increasing toward the bottom of the panel. The inset in each case zooms out to reveal the height of the lineshapes 
for all the momenta. While the linewidths near kp are strongly effected by rising temperature, the incoherent parts particularly 
at high frequency have very little temperature dependence. Note also that the tails in Case (A) demonstrate an inte rest ing 
structure near uj — —AeV, a second maximum far from the nominal QP. A similar structure was previously observed irP^. It 
appears in Case (A) but not in Case (B). Case (B) displays a local minimum of the linewidths for /c > /cf, as seen more clearly 
in Fig. (IT7). 



a parameter A to second order yields the equations studied here, and has been argued to be the first step in a 
systematic microscopic study of the t-J model. The range of validity of this calculation is somewhat removed from 
the most interesting regime of optimal doping, and is self consistently determined to be n<0.7. This theory differs 
from the simplified ECFL theory, used to model the optimal doping system line shapes in Ref. (2), in two respects. 
In the present work the various self energies are calculated self consistently starting from the t-J model rather 
than postulated- indeed the assumed forms of the self energies turn out to be quite close to the present calculations, 
thereby validating the approximation. Secondly, the static multiplicative renormalization of the bare band energies 
Ck in Eq. ([5| is also somewhat different from that in Ref. (2); it is chosen here by a systematic expansion in A as 
formulated in Ref. (j5]). We must therefore be careful in comparing the present results with data at higher particle 
densities, e.g. closer to optimal doping. A higher order in A calculation, than the O(A^) calculation presented here, is 
needed for addressing higher particle densities. Nevertheless, our work represents the first systematic and essentially 
analytical calculation of the t-J model at the above mentioned densities, and is interesting for that reason. It is 
encouraging to note the broad similarity of our results with other numerical works mentioned in the introduction and 
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FIG. 17: For Cases (A) and (B) at (n, T) = (.75, 130K) the decay rate calculated from the Dysonian self energy: F^ = pi:{k, Ek) 
is plotted in the first quadrant of the Brillioun zone. We see a second set of minima, particularly prominent in Case (B), apart 
from the expected minimum at the FS. Note also that the value of F at the FS is much lower in Case (A) than in Case (B). 
This occurs because in Case (B) the chemical potential sits near a van Hove singularity resulting in high scattering. 
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FIG. 18: For {n,T) = (.75, 130K) the spectral function pg is plotted for k's along the line perpendicular to the anti-nodal FS. 
For Case (A) on the left and Case (B) on the right this corresponds to the (01) direction and the line from (0, tt) (vr,7r), 
respectively. The bandwidths seen along these lines are much smaller than those seen along along the nodal line. Thus the 
momentum dependence of the self energies in the present theory distinguish significantly between the different directions. 



text above. 

In Section (III) we have presented detailed results for the chemical potential and the momentum distribution 
function m/c, for two standard sets of band parameters of current interest. We have then detailed the various energy 
dispersions that emerge in our calculation- these include the bare band dispersion, the MDC dispersion and the EDO 
dispersion. We also evaluated some measures of the line asymmetry. The so called kink feature is displayed and 
its dependence on the band parameters is explored. The wave function renormalization is evaluated at various 
densities and the Fermi velocity computed along various directions. We observed a second (shadow Fermi) surface, 
outside the usual Fermi surface, representing a shallow minimum in the quasiparticle relaxation rate. It is different 
from the usual Fermi surface in location, size and also in the non zero (but small) dissipation at the surface. It is 
more pronounced for Case (B) than in Case (A), but not very visible in the momentum space distribution function 
rrik in either case. 

The relaxation rates of quasiparticles at finite T and resistivity as well as the local DOS are computed. We displayed 
the spectral function at various densities and temperatures along different directions in detail. We also explored the 
effect of elastic impurity scattering on the line shapes in detail. Finally we displayed the optical conductivity and 
resistivity from this theory, after throwing out vertex corrections, and compared with some data on overdoped systems 
Ref. (25). As far as possible, all computed objects are presented on an absolute scale^ using lattice parameters typical 
of the Cuprate materials, and a fundamental energy scale set by the nearest neighbor hopping t ~ 3000K. It is clear 
that one could tune this scale to better fit a given experiment, e.g. the optical conductivity in Fig. (21) gives a good 
fit on an absolute scale to Case (A) if we rescale t, whereas it does well in Case (B) with the present scale. 
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FIG. 19: (n, T) = (.75, 130K) The effect of the elastic scattering term 77 on the spectral functions pg along < 11 > in Cases (A) 
and (B) (left and right, respectively). For comparison we also include the pristine system with 77 ^ 0. and display the resulting 
curves for three values ?7 = 0, .032, .12 eV. The two nonzero values are typical for laser and synchrotron ARPES respectively, 
as argued in Ref. ((2|. A large (synchrotron ARPES) value of 77 ^ .12 eV leads to extended tails and a greatly reduced peak 
value of the spectra, while the smaller (laser ARPES) value of 77 ^ .032 has a less drastic effect on the pristine case. 
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FIG. 20: Case (A) n = .75, T = 60,90, 130, 190, 280, 410, 605i^^. The optical conductivity is calculated (on an absolute scale) 
and illustrates how increasing T "fills up" the regime < uj < 1000 cm~^ rapidly. The rise of conductivity at very low cj is 
also inferred from the DC resistivity displayed in Fig. |l3t. The phase of the complex a falls off rapidly beyond 4000 cm~^. 
Results for Case (B) are similar and shown below in Fig. (|21|). 



R. R. P. Singh, R. L. Glenister, Phys. Rev. B 46, 14313 (1992). 
M. M. Zemljic and P. Prelovsek, Phys. Rev. B 72, 075108 (2005). 

M. M. Zemljic, P. Prelovsek, and T. Tohyama, Phys. Rev. B 76, 012502 (2007), M. M. Zemljic, P. Prelovsek, and T. 

Tohyama, Phys. Rev. Letts. 100, 036402 (2008). 
^2 J. Kokalj and P. Prelovsek, Phys. Rev. B 75, 045111 (2007). 

J. Jaklic and P. Prelovsek, Adv. Phys. 49, 1 (2000). 

C. Gros, R. Joynt, T.M. Rice, Phys. Rev. B 36, 381 (1986). 

B. Edegger, V. N. Muthukumar and C. Gros, Adv. Phys. 56 927 (2007). 
^'^ A. Paramekanti, M. Randeria, and N. Trivedi, Phys. Rev. Letts. 87, 217002 ( 2001). 

B. Moritz, F. Schmitt, W. Meevasana, S. Johnston, E. M. Motoyama, M. Greven, D. H. Lu, C. Kim, R. T. Scalettar, Z.-X. 

Shen, T. P. Devereaux, New J. Phys. 11, 093020 (2009). 

W. Metzner and D. Vollhardt, Phys. Rev. Letts. 62, 324 (1989). 
2° A. Georges, G. Kothar, W. Krauth, and M. J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 

X. Deng, J. Mravlje, R. Zitko, M. Ferrero, G. Kotliar, A. Georges, arXiv:1210.1769 (2012). 

W. H. Press, S. A. Teukolsky,W. T. Vetterling, B. P. Flannery, Numerical Recipes: The Art of Scientific Computing, 
Cambridge University Press, New York, NY (1988) 
2^ T. Valla, A. V. Fedorov, P. D. Johnson, Q. Li, G. D. Gu, N. Koshizuka, Phys. Rev. Letts. 85, 828 (2000) 
J. Graf, G.-H. Gweon, and A. Lanzara, Physica C: Superconductivity, 460-462, (2007) 
A. V. Puchkov, D. Basov and T. Timusk, J. Phys.: Condens. Matter 8, 10049 (1996) 



21 




Case (A) 
r=130,190,280ir 



500 1000 1500 2000 

(j) (cm"^) 




500 1000 1500 2000 

(jL> (cm"^) 




2000 



FIG. 21: We make an explicit comparison here with optical conductivity measurement in Puchkov et al in Ref. (|25p . The 
experimental data is reproduced from Ref. (|25|) with the authors kind permission. It pertains to an overdoped Thallium based 
compound with Tc = with a density that is expected to be close to the value assumed in the calculation n ^ .75. We note 
the similarity of magnitude and scale with our Case (B), our Case (A) leads to a bigger value of Ylea{uj). 



